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The study of proton-antiproton interactions with an 
isolated, high energy photon and an identified 6-quark jet 
is a testing ground for quantum chromodynamics (QCD) 
predictions at the Tevatron. At photon transverse en- 
ergies Ej, below 70 GeV Compton scattering processes 
dominate production, with gb — > *ygb or qb — > jqb domi- 
nating at lower Ej depending on the relative sizes of the 
quark and gluon parton density functions. Above an Ej, 
of 70 GeV qq — > bbj quark annhiliation processes domi- 
nate production A cross section measurement there- 
fore provides a well measured probe of the hard scatter- 
ing dynamics within the proton, and some sensitivity to 
the 6-quark content of the proton, whose parton density 
function is so far indirectly extracted from constraints on 
the gluon density functions. 

The first measurement of photon and heavy flavour jet 
production was performed using 86 pb _1 of integrated 
luminosity taken at ^/s = 1.8 TeV with the CDF I de- 
tector 0. Heavy jet flavor was signalled by muons con- 
tained in the jet. The results were interpreted in terms 
of new physics processes involving the radiative decay 
of Techni-Omega states [H, and radiative decays of su- 
persymmetric particles in gauge-mediated supersymmc- 
try models [4|. Recently the DO collaboration has mea- 
sured the cross section of heavy flavour jet production 
in association with a photon [5|] using data collected at 
=1.96 TeV. In this paper we exploit improvements in 
the CDF II detector to identify b jets by a lifetime based 
secondary vertex tag, use a larger dataset collected at a 
higher center-of-mass energy, probe lower photon trans- 
verse energies, and employ a superior analysis technique 
where all backgrounds are determined directly from data. 
We present our results as the differential photon + b jet 
production cross section, as a function of photon Ej, and 
the total photon + b jet production cross section. 

The CDF II detector is described in detail in [6j . Only 
the components which are most relevant to this analysis 
will be described here. 

The detector is composed of a central spectrometer in- 
side a 1.4T magnetic field, surrounded by electromagnetic 
and hadronic calorimetry and muon chambers. The spec- 
trometer is composed of a multi-layer silicon vertex de- 
tector inside a cylindrical multi-wire drift chamber. The 
combination of these tracking detectors measures charged 
particle trajectories with a transverse momentum (pt) 
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precision of Apx/pj = 0.07% (GeV/c) -1 , and an uncer- 
tainty on the transverse impact parameter of about 40 
fim for tracks of pr above 1 GeV/c, which includes the 
intrinisic beam size of about 30 /im. Information from 
the central tracker can be sent within trigger latency to 
the hardware tracker SVTpj], that provides an extra- fast 
measurement of the tracking parameters the can be used 
at trigger level. 

Segmented sampling calorimeters, arranged in a pro- 
jective tower geometry, surround the tracking detec- 
tors and provide energy measurements within the region 
| r\ |< 3.6. Central calorimeters [8[ cover the region 
| rj |< 1.1, with an electromagnetic (hadronic) energy 
resolution of a(E)/E = 13.5%/ Vl? © 2.0% {a(E)/E = 
50% I \/Et © 3%). The end- wall hadronic calorimeter ex- 
tends this coverage to | r\ |< 1.3 @ with an energy resolu- 
tion of 75%/ v / £r © 4%, whilst the region 1.3 <| 77 |< 3.6 
(1.1 <| 77 |< 3.6 for the electromagnetic calorimeter) 
is covered by forward calorimeters [10j , with hadronic 
and electromagnetic energy resolutions of 80% /V~E © 5% 
and 16%/ \f~E~ © 1% respectively. In order to distinguish 
electromagnetic clusters from photons and electrons, and 
clusters from neutral pion decays, the central electromag- 
netic calorimeter is equipped with two sets of wire cham- 
bers: a preshower detector (CPR) situated in front of 
the active part of the calorimeter to detect early photon 
conversions in the solenoid coil, and a shower maximum 
detector (CES) placed at the position of maximal width 
of electromagnetic showers to measure the shower profile. 
Signals from these detectors are compared with expected 
deposits from electromagnetic clusters, obtained from 
testbeam data. Each electromagnetic cluster is given a 
weight related to the probability that it originates from 
a photon. 

We use data obtained by two triggers; a trigger which 
requires a photon-like object with transverse energy 
larger than the threshold for the inclusive photon trig- 
ger at CDF of 25 GeV ('high Ej photon'), and a trigger 
('SVT photon') which requires a photon-like object with 
transverse energy larger than 12 GeV, a jet with trans- 
verse energy larger than 10 GeV, and a track, measured 
by the SVT Q, with transverse momentum larger than 
2 GeV/c, an absolute rapidity smaller than 1 and an im- 
pact parameter larger than 120 /im. The SVT compares 
hits from the tracking detectors with pre-fitted tracks 
stored in an associative memory to extract track param- 
eters in real time. In this way tracks with significant im- 
pact parameters can be quickly identified. The impact 
parameter resolution of this procedure is less than 50 pui 
(including the contribution from the beam size), similar 
to the accuracy obtained with full detector reconstruc- 
tion. An integrated luminosity of 340 (208) pb~ x of data 
has been analyzed in the high Ej, photon (SVT photon) 
triggered dataset. 

The two triggers have very different efficiencies. The 
high Ej photon trigger has an efficiency close to 100% 
for events with Ej above 28 GeV. To obtain photons of 
lower Ej, the SVT trigger must be used. The efficiency 
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of the SVT photon trigger is found by comparing the 
number of events that pass both triggers. Being based 
on strong requirements on a single track rather than a 
more refined secondary vertex reconstruction, the effi- 
ciency of the SVT-based photon trigger is about 50% 
(with a relative uncertainty of 4%, due to the size of the 
overlap sample) and was found to be independent of the 
jet transverse energy. 

Selected candidate events must pass one of the two 
photon triggers, contain an isolated central (| r\ |< 1.1) 
photon of E T > 20 GeV, and a b jet of E T > 20 GeV 
within | r) |< 1.5. 

Photon candidates must have a calorimeter cluster 
which contains predominately electromagnetic energy. 
The cluster hadronic energy fraction Eh ad / E^m must 
satisfy E HAD /E EM < 0.055 + 0.00045 * ET>, where is 
the photon energy and Eh ad and Eem are the hadronic 
and electromagnetic energy deposits within the cluster. 
The electromagnetic shower profile must also agree with 
that expected for an electromagnetic deposit. In or- 
der to reduce contamination from neutral meson (such 
as 7T°) decays, photon candidates must be isolated from 
nearby calorimetric deposits and tracks. We require that 
E T (R 0A ) < 2.0 + 0.02(£j - 20), where E T (R 0A ) is 
the summed calorimetric cluster transverse energy de- 
posits in a cone of radius R = y A0 2 + Arj 2 = 0.4 
around the photon candidate, to ensure calorimetric iso- 
lation. In addition, the summed track transverse mo- 
menta of tracks inside this cone, pt(R 0A ) must satisfy 
p T (R 0A ) < 2.0 + E T * 0.005 to ensure isolation in the 
tracking detectors. Events containing adjacent calori- 
metric clusters in the CES, which can signal the presence 
of neutral mesons, are rejected. 

The jets in each event are reconstructed using the JET- 
CLU algorithm [U with cone radius R = 0.4 (0.7) for 
events containing photons of E T < (>)26 GeV. To re- 
cover the true parton energy, jets are corrected for in- 
strumental effects 12j. We select events containing at 
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least one jet with corrected Et > 20 GeV, and whose 
axis lies outside a cone of AR = 0.7 surrounding the 
photon candidate. In order to identify jets arising from 
b hadrons we search for displaced secondary vertices [l3| . 
A b jet is identified ('tagged') when the secondary ver- 
tex is more than two standard deviations away from the 
beam position, and in the same direction away from the 
beam position as the jet momentum. We require at least 
one b jet to be identified in each event. The efficiency 
of the ^-tagging algorithm is just under 25% for b jets of 
E T = 20 GeV, and increases to 40% at E T = 50 GeV. 

The PYTHIA [3] Monte Carlo code is used to simulate 
photon + jet production at leading order and to estimate 
the photon and jet selection efficiencies. The Q 2 scale of 
the interaction is set to 225 GeV 2 , and the CTEQ5L [l5[ 
parton distribution functions are used. A simulation of 
the underlying event is included (l6| . 

Backgrounds to photons arise from high energy ir s, 
which decay to pairs of overlapping photons that can- 
not be distinguished in the calorimeter, but whose rate 
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FIG. 1: Fit to the invariant mass of tracks composing the sec- 
ondary vertex in data, for photon candidates having Et > 26 
GeV. The points are data, and the stacked, shaded histograms 
represent the estimated contributions of the b, c- and light- 
quark jets. 



can be estimated directly from the data. The fraction 
of misidentified photons is estimated using the clusters 
detected by the CPR and CES detectors of the electro- 
magnetic calorimeter. Photon pairs produced by 7r° de- 
cays will produce a larger number of hits in these de- 
tectors, and each photon candidate is assigned a weight 
proportional to the probability of being a signal or back- 
ground cluster following the procedure detailed in [l7l |. 
The fraction of correctly identified photons in the sam- 
ple passing those selection criteria increases with increas- 
ing Et, going from about 50% at the lower end of the 
spectrum considered here (20 GeV) to around 80% at 
high E T . Backgrounds to b jets can arise from jets orig- 
inating from charm quarks, since charm hadrons have 
a lifetime between a quarter and two thirds that of b- 
hadrons, and jets originating from light quarks and glu- 
ons, where random combinations of tracks can sometimes 
mimic a vertex displaced from the beamline. The purity 
of our selected events is determined by fitting the invari- 
ant mass of tracks composing the secondary vertex in 
data to Monte Carlo templates of the shapes expected for 
b-, charm (c-) and light-quark jets. Fig. [T] shows an ex- 
ample of the fit to data. Here, about a third of jets arise 
from b quarks. This invariant mass is generally lower 
than the corresponding hadron mass due to mis-assigned 
tracks and unreconstructed neutral hadrons. It can be 
seen that the template shapes of the different quark jet 
types are sufficiently different to provide reasonable dis- 
criminating power. 

Identified b jets can arise from either the photon + 
b jet signal, or from events where one or both objects 
(photon and jet) are mis-identified. We assume that the 
tagged jet composition in ir° + tagged jet events is similar 
to 7r ± + tagged jet events. We estimate the b jet purity of 
the misidentified photon background using data collected 
with a jet based trigger. Events are required to contain 
two jets, one of which must be tagged and have sim- 
ilar transverse energy and pseudorapidity requirements 
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to the b jet, and a second which passes similar kinematic 
requirements to the photon. The fraction of b jets in 
this sample can be found by fitting the secondary vertex 
invariant mass of the tagged jets. The purity of the se- 
lected jets ranges from 50% for jets of Et around 20 GeV, 
to about 15% for jets of Et around 75 GeV where the 
rate of light-quark jet tagging increases. This b fraction 
is then normalised to the estimated number of misidenti- 
fied photons, and subtracted from the estimated number 
of b jets in the whole event sample to yield the total 
number of photon + b jet events. 

Some 10900 (55800) events pass the selection criteria 
in the high transverse energy photon (SVT photon) trig- 
gered datasets. Candidate events are divided into bins of 
photon transverse energy. The numbers of events in each 
bin are corrected for background contributions, trigger, 
selection and acceptance efficiency, and divided by the 
appropriate integrated luminosity to yield a cross sec- 
tion. Results are given in Table U which also lists the 
systematic uncertainties detailed later. Note that the 
statistical uncertainty for the high E2, photon dataset 
includes contributions from finite Monte Carlo statistics. 

Several sources of systematic uncertainty in the cross 
section determination have been studied: photon iden- 
tification, jet energy scale, b jet identification, and lu- 
minosity. In the following only the largest contributions 
will be quantified. 

Variables used in photon identification have been stud- 
ied in Z — > e + e~ simulation and data events [l8|. Sim- 
ulation and data are consistent. Uncertainties in the 
misidentified photon estimate arise from assumed values 
for the preshower detector hit rate, the rate of backscat- 
tered showers and the fractional composition of fake pho- 
ton backgrounds. The associated systematic uncertainty 
is about 6%. 

The systematic uncertainty on the jet energy scale has 
been studied in detail elsewhere [l2j and the findings ap- 
plied to this analysis. It decreases with increasing jet 
Et, and is about 5% for jets of 35 GeV Et- Uncertain- 
ties have also been determined for the effect of multiple 
interactions overlapping in the same data event, and the 
additional uncertainty on the jet scale introduced by a b 
jet. 

Uncertainties in the b quark purity of the sample arise 
from imperfect modelling of the template shapes used in 
fits to extract the b purity. Differences in shape can arise 
if the secondary vertex invariant mass of jets containing 
one ('single') or two ('double') b quarks differs, or if track 
efficiency is incorrectly modelled. We find that template 
shapes obtained from single and double quark jets are not 
consistent with each other. The data are fitted to tem- 
plates composed of mixtures of the individual single and 
double quark templates (ranging from to 100%), and 
the x 2 of the resulting distributions with respect to the 
data is computed. We take as a 1 — a deviation the value 
for this mixture for which the \ 2 increases by one unit 
with respect to its minimum value, and recalculate the 
cross section using this mixed template. The systematic 



uncertainty is the difference between this value and the 
default cross section (obtained with the default diquark 
fraction) . 

The uncertainty on the cross section is assigned as the 
difference between the reference value and that obtained 
from the value of the single/double quark template for 
which the \ 2 increases by one unit with respect to its 
minimum value. This is the largest single source of un- 
certainty and is about 17%. Previous studies [n| sug- 
gest a difference in tracking efficiency between data and 
simulation which is a function of isolation, momentum, 
and position in the detector. We remake the invari- 
ant mass templates incorporating the inefficiency derived 
from data, and take the full difference between results ob- 
tained using these and the nominal templates (5%) as the 
systematic uncertainty due to this source. 

Other systematic uncertainties on b jet identification 
arise from the difference in measured tagging efficiency 
between data and simulation, from the difference in the 
tagging efficiency of single and double b jets, and from 
simulated b hadron multiplicity. The first uncertainty 
has been previously studied [l3| ■ The difference in scale 
between tagging efficiency in data and simulation was 
found to be 0.91 ± 0.06. This results in a 6% uncertainty 
on the measured cross-sections. The uncertainty on tag- 
ging efficiency for single and double b jets is determined 
as the difference between results obtained using the frac- 
tions of single and double quark templates corresponding 
to ±1 standard deviation, as found earlier, and is about 
7%. We have adopted the findings of previous studies 
of the effect of assumed b hadron multiplicity [l9[ (a 1% 
effect on the measured cross section). The SVT-based 
analysis is also affected by the statistical precision of the 
trigger efficiency determination (about 10%). Finally, the 
luminosity is subject to a ±6% uncertainty (20| . 

The cross section for photons produced in association 
with b jets is tabulated in Table fl] Both differential 
and inclusive results are given, with their statistical and 
systematic uncertainties. Note that results from both 
datasets are listed. There is overlap between the final 
SVT photon and the first high E T photon bins, which 
give consistent results. Due to its greater statistical pre- 
cision, the former bin is used in the final results. The 
measurements are corrected to hadron level [21| , so that 
they can be directly compared to a recent next-to-leading 
order calculation [lj. This prediction was derived analyti- 
cally, using the CTEQ6.6M parton density functions [22j |. 
and a renormalisation, factorisation and fragmentation 
scale set to the transverse momentum of the photon. It 
does not include non-perturbative effects (hadronisation 
and underlying event), and is presented in terms of par- 
ton level jets. 

The measured cross sections are shown compared with 
this prediction in Fig. [5] Also shown are the theoreti- 
cal uncertainties due to (i) choice of scale and (ii) uncer- 
tainty in parton density function assumption. Agreement 
with next-to-leading order theory is good over the entire 
photon E T range probed. It should be noted that due 
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E T (GeV) da(pp^ 7+ > lb-jet)/rfE%(Data) d<r(pp 7+ > lfe-jet)/rfgj.(NLO) 



20 


- 24 


3 90 + 49 + 84 


3 97 + 78 


24 


- 28 


3.01 ± 0.41 ± 0.63 


3.67 ± 0.32 


26 - 


28(*) 


3.13 ±0.51 ±0.67 


3.01 ±0.21 


28 


- 31 


2.90 ±0.42 ±0.61 


2.65 ±0.18 


31 


- 35 


1.24 ±0.20 ±0.27 


1.72 ±0.14 


35 


- 43 


0.94 ± 0.14lg;2o 


0.92 ±0.10 


43 


- 70 


0.20 ± 0.03 ± 0.04 


0.21 ±0.05 



TABLE I: The measured differential cross section for central photon production in association with at least one &-jet of Et > 20 
GeV, inside | r\ j< 1.5, tabulated as a function of the photon transverse energy E T . The first (second) uncertainty quoted is 
statistical (systematic). Results are compared to a next-to-leading order (NLO) prediction, whose uncertainty arises from the 
numerical integration procedure [ljj. Note that the first two measurements are made using the SVT dataset, and the remainder 
with the high E T dataset. The overlap bin, denoted by (*), is not used in the final results. Systematics about luminosity, 
secondary vertex tagging scale factor, the effect of multiple verteces are correlated between all bins. The uncertainty due to 
the statistical precision of the SVT trigger efficiency is fully correlated between the first two bins only. All other uncertainties 
are uncorrelated between bins. 




Photon ^ (GeV) 



FIG. 2: Top: b ± photon cross section as a function of photon 
Et, compared to NLO QCD calculations. The light dashing 
is a quadrature sum of uncertainties coming from scale varia- 
tion (both renormalisation and factorisation scales varied by 
a factor 2 and 0.5) and PDFs, while the darker dashing rep- 
resents the scale variation contribution only. The inner error 
bars for data represent the statistical uncertaintiess, the outer 
the combination of statistical and sytematic. The bottom plot 
shows the ratio of data to NLO calculation, where the error 
bars and shading have the same meaning as before. 

to numerical stability problems, the first bin in the NLO 
calculation starts at 18 GeV instead of 20 as for the data. 

The total cross section <r(pp — > 7+ > 16-jet; E T > 20 
GeV) has been measured to be 54.22 ± 3.26 (stat) ± 5 5 ™ 
(syst) pb. This is consistent with the next-to-leading 
order prediction of 55.62 ± 3.87 pb. 

In conclusion, the cross section for photon production 



in association with b jets has been measured in proton 
antiproton collisions at y/s — 1.96 TeV with the CDF 
II detector. The measurement has been made for &-jets 
with Et > 20 GeV inside | 77 |< 1.5, and for photons 
of at least E T > 20 GeV inside | 77 |< 1.1, including the 
lowest photon transverse energies probed to date. The re- 
sults are consistent with next-to-leading order theoretical 
peturbative QCD predictions, using CTEQ6.6M parton 
density functions, throughout the photon E T range mea- 
sured, while leading-order calculations would predict a 
cross section smaller by about 30%. The level of accu- 
racy of this measurement is therefore already sufficient 
to discriminate between the first orders of perturbative 
expansion, and favor the most precise NLO predictions. 
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